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ABSTRACT 

Optimum  locations  were  found  for  dovap  and  radar-ranging 
systems.  The  locations  of  ground  stations  given  are  those  which, 
for  a fixed  missile  position,  minimize  the  effect  of  the  geometry 
of  the  configuration  on  the  precision  of  the  determination  of 
missile  position.  Certain  optimum  locations  arc  also  given  for 
hyperdop,  which  are  not  usable  on  a ballistics  test  range. 
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INTRODUCTION 

Optimum  locations  were  found  for  dovap  and  radar-ranging  systems  in  the 
investigation  reported  here,  extending  considerably  the  work  done  in  Ref.  1 and 
2.  The  investigation  was  concerned  solelv  with  the  mathematical  problem  as 
defined  in  these  references,  The  locations  of  ground  installations  given  are 
those  which,  for  a fixed  missile  position,  minimize  the  effect  of  the  geometry 
of  the  configuration  on  the  precision  of  the  determination  of  missile  position. 

Basic  results  and  background  of  the  optimum-location  problem  are  given  in 
Ref.  1,  2,  and  3 and  are,  therefore,  used  without  discussion.  As  in  the  references, 
the  volume  of  the  ellipsoid  of  concentration  is  used  in  this  report  as  a measure 
of  precision.  The  results  of  Ref.  1 and  2 that  the  geometry  of  the  configuration 
has  an  effect  on  precision  quite  separate  from  other  factors  are  also  used.  It  is 
evident  from  these  results  that  only  the  angle’  between  the  several  rays  from 
ground  installations  to  missile  position  are  involved  in  the  effect  of  the  geome- 
try, when  the  assumptions  given  and  discussed  in  Ref.  1 and  2 are  satisfied  and 
when  the  stations  are  given  equal  weight  factors  (P;  or  Pj;-).  These  assumptions 
are  made  and  these  conclusions  used  throughout  this  report.  The  volume  of  the 
ellipsoid  of  concentration  is  proportional  to  o3/\/S,  where  o is  the  common 
standard  deviation  of  the  observations  and  where  A is  a certain  determinant 
determined  by  these  angles.  Ar.  optimum  location,  then,  is  one  which  yields  an 
absolute  maximum  of  A for  a fixed  number  of  ground  stations. 

The  expressions  given  in  Ret.  1 and  2 for  the  determinant  for  each  of  the 
three  types  of  instrumentation  are  easily  seen  to  involve  only  the  direction  co- 
sines of  the  rays  from  ground  station  to  “approximate  missile  position”  and  not 
the  slant  ranges  from  station  to  this  position.  Since  the  sum  of  squares  of  the 
members  of  a set  of  direction  cooines  is  always  exactly  unity,  A is  always 
bounded  and  has  a maximum.  It  is  convenient  in  the  present  report  to  arrange 
these  direction  cosines  in  the  3-row  matrix  X ■=  (x,-,),  in  which  the  /th  column  is 
the  set  of  direction  cosines  of  the  ray  from  the  /th  ground  station  to  the  missile, 
for  / ■»  J,  2,  • • • , n.  For  each  of  the  types  of  instrumentation,  it  is  shown  below 
that  there  is  a matrix  A of  constants,  determined  by  the  number  of  stations  and 
the  type  of  instrumentation,  and  such  that  A.  = \XaX'\,  where  A"  is  the  transpose 
of  A’. 

The  results  obtained  in  the  present  investigation  are  first  summarized  and 
then  derived  in  the  separate  sections  below  devoted  to  the  respective  types  of 
instrumentation.  The  treatment  of  the  radar  case  follows  the  treatment  of  Ref.  2 
in  not  making  use  of  azimuth  and  elevation  information.  (When  Ref.  2 was  pre- 
pared, the  equipment  available  did  not  vield  adequate  information  of  this  kind.) 
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Optimum  locations  for  any  number  of  stations  are  given  for  radar  and  dovap. 
Certain  locations  are  shown  to  be  singular  for  hyperdop  in  that  the  data-reduc'.ion 
methods  of  Ref.  1 fail,  while  the  only  optimum  locations  found  for  hyperdop  are 
not  usable  on  a ballistics  test  range.  Numerical  results  showing  the  effect  on 
accuracy  of  the  number  of  stations  and  f.he  vrriation  in  accuracy  as  the  missile 
follows  a trajectory  through  the  optimum  are  given  for  radar  and  dovap.  In  con- 
clusion, the  report  offers  some  general  discussions  on  certain  results  found. 

Appendix  A gives  mathematical  proofs  of  certain  theorems  relating  to  deter- 
minants of  the  form  \XAX'\.  Two  of  these  theorems  are  essential  in  the  treat- 
ments of  the  three  types  of  instrumentation: 

THEOREM  1.  Rigid  rotations  in  spacer  of  any  configuration  leave  A unchanged. 

ItaEORlM  2.  The  determinant  A is  never  greater  than  the  product  of  its  main 
diagonal  elements  and  is  equal  to  this  product  if  and  only  if  the  latter  is  zero  or 
all  the  off-diagonal  elements  are  zero. 

THE  RADA?  CASE 


Summary 

The  known  optimum  locations  in  the  radar  case  may  be  arranged  for  con- 
venience in  two  series,  both  of  which  involve  circular  cones  with  the  missile 
at  the  vertex.  Because  of  Theorem  1,  these  cones  may  be  considered  initially 
to  have  vertical  axes. 

Tbe  first  and  more  important  series  has  a fixed  element-axis  angle  of  arc 
cos  1/V3~54°44'  and  the  stations  are  on  elements  of  the  cone.  A case  for 
five  stations  is  illustrated  in  Fig.  1.  The  station-missile  rays  may  be  (1)  equally 
spaced  around  the  cone  or  (2)  divided  into  groups  of  not  less  than  three  stations 
each  with  the  static. is  of  a group  equally  spaced  and  tbe  groups  arbitrarily 
oriented  with  respect  to  each  other.  These  are  -he  only  spacings  in  this  series 
for  n < 6 stations:  other  ■iwrino*  «r*  possible,  hctvsvcr,  lot  larger  Values  of 
n (shown  in  Appendix  B), 

The  second  series  has  some  stations  on  the  axis  of  the  cone.  A construc- 
tion in  this  srrics  is  possible  only  if  the  nuud.-r.  *,  of  stations  on  tbe  axis  is 
not  greater  than  nf 3,  where  n is  the  total  number  of  stations.  If  r » n — s,  the 
element-axis  angle  of  the  cone  is  arc  cos  y'(l/r)[(n/3) sj.  The  r stations  not 
on  the  axis  are  on  elements  of  the  cone,  spaced  about  it  according  to  the  same 
rules  as  used  for  the  first  series. 

In  addition  to  Theorem  1,  another  rule  for  deriving  further  optimum  locations 
from  a known  optimum  location  is  available  in  the  radar  case.  The  rule  is  that 
changing  the  sense  of  a station-missile  ray  leaves  A unchanged. 

These  results  are  derived  in  detail  below. 

Derivation 

The  first  result  proved  in  this  section  is  the  rule  stated  above  for  deriving 
further  optimum  locations.  This  result  is  easily  derived  from  the  expressions  for 


2 


NAVORD  REPORT  2099 


RANGE  PLANE 

FIG.  I.  Typical  Configuration  of  Stations  said  Missile. 


the  elements  of  A.  Such  expressions  are  given  on  page  3 of  Kef.  2 in  terms  of 
quantities  (Xj,  and  yj,  which  are  the  direction  cosines  x.;,  r;;s  and  of  the 
present  report.  Thus,  expression  9 on  page  3 of  Kef.  2 is  simply  X XX' (A  is 
the  unit  matrix),  and  the  element  in  the  CXth  row  and  Bih  column  of  A is 

n 

Z *cu*/3> 

t"l 

Evidently,  this  element  is  unchanged  if  xj,-,  x2l-,  and  *3,-  are  replaced  by  -Xu, 
— x2l-,  and  —13,',  respectively;  hence,  the  rule  is  obtained. 

The  second  result  proved  is  that  locations  of  the  first  series  are  optimum. 
To  show  tnis,  note  that,  when  the  cone  axis  is  vertical,  the  vertical  components 
x2,  are  al!  equal  and  have  the  value  ]/\/3. 

With  the  stations  in  the  range  plane  and  the  station-missile  rays  equally 
spaced  around  the  cone,  the  stations  are  then  equally  spaced  in  the  range  plane 
around  a circle  whose  center  is  directly  below  the  missile.  The  pairs  (xj;,  X3;) 
are  direction  numbers  in  the  range  plane  of  the  projections  into  this  plane  of 
the  station-missile  rays.  Since  X\t  + x|;  «=  1 - x2l-  = 2/3,  x'u  = 3/2  and 
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*3i  = *3 iV^/2  are  direction  cosines  in  the  range  plane  of  these  projections.  The 
complex  numbers  z;  = x\i  + X3,V-1  lie  on  the  unit  circle  in  the  complex  plane 
and  are  equally  spaced  about  this  circle;  hence,  they  are  the  nth  roots  of  some 
complex  number  c<p  whose  absi  lute  value  is  unity.  Thus,  zj,  • • - , z„ 

are  the  roots  of  the  equation 

(1)  z"  -exp  (0-/-1;  - 0 (n  > 3) 

However,  it  is  well  known  for  any  polynomial  of  the  form 
Cn+alCn~1  + a,Cn-2^‘“  + an 
with  zeros  £]>  •••,£„  *tat  the  relations 

O'  Zb"* 

I te,  - «! 


hold.  From  the  first  of  these,  since  a;  = 0 in  the  polynomial  in  Eq.l, 

Z (*li  + *31 V^)  « 0 or  S*li  = 0 und  £ *3i  = 0 

and  it  follows  easily  thut 

(3)  £*ii*2£  = (l/v/3)  £ Xu-  = 0 and  £ *3«*2i  = <l/\/3)  £ *3i  = 0 

From  the  second  of  relations  2,  since  a2  = 0 in  Ea.  1, 


«*/  = <i/2)(  Z Z *,*,  -±>?)  — n/2)  ±zf  = 0 

\ i-1  /-I  /-!  / 1-1 


where  the  second  step  follows  from  X z;-  = 0.  However, 


£ */  = £ t(*iy)2  - (*3/)2l  + 2\Z-i  £ *i/*3/ 
;*i  y-i  ;-i 


so  that 


X *l/*3/  = 0 


Z *?,  = Z *3/ 


Relations  3 and  4 show  that  the  off-diagonal  elements  of  |/VA’,|  vanish  so 
that  A is  equal  to  the  product  of  its  main  diagonal  elements.  If  these  are 


Scx=  £ *h 

i-l 


fa  = 1,  2,  3) 


then,  using  x\\  + x2i  + x32,-  = 1, 

(6)  iq  + S2  4-  S3  = n 

and,  since  x2;  - 1/V3,  S2  =>  n/3  and,  from  Eq.5,  Sj  = S3  al3o.  This  means  tout 
A = S1S2S3  = (n/3)3. On  the  otherhand,  the  partial  derivatives  of  Q = S2(n  - A'i  — b2) 


dQ 

— = S2(n  - 26’i  - S2) 
oS; 
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dQ 

= St(n  - Sj  - 2%) 

Caj2 

The  obvious  nontrivia!  solution  of  the  equations  setting  these  derivatives  equal 
to  0 is  Sj  = S'2  «•  n/3.  Thus,  (n/3)3  is  the  absolute  maximum  of  the  product 
which  A never  exceeds  by  Theorem  2.  The  configurations  of  the  first  series  for 
which  the  stations  are  equally  spaced  around  the  cone  thus  yield  the  maximum 
value  of  A and,  hence,  are  optimum. 


For  the  configurations  of  the  first  serieo  in  which  the  stations  are  grouped 
into  r group.i,  the  argument  can  be  repeated  in  its  entirety  except  that  the  poly- 
nomial of  Eq.  1 becomes 


[z"»  - exp  - exp  (fljv'-l)]  • • • lz”'-  exp  ((9rV/-l)] 

where  each  n*  > 3 and  a\  = 02  = ® 38  before.  The  other  locations  not  specified 
above  fellow  from  further  results  involving  complex  polynomials  given  in 
Append :x  B. 


The  final  result  proved  in  this  section  is  that  locations  of  tit*'  second  series 
are  optimum.  In  this  case,  the  matrix  A takes  the  form 

I *hVl  - y2  *i2\A  - y2  • • • x'u'A  - y2  o . . . o\ 

I y y •••  y 1 • • • 1 I 

\ *3iVl  - y2  *32\/l  - y2  • • • - y2  0 • • • 0 / 

where  y is  the  elevation  component  common  to  all  stations  on  the  cone,  where 
the  last  s columns  of  X apply  to  the  stations  on  the  cone  axis,  and  where  xjj  and 
r3J  may  he  taken  as  the  real  and  imaginary  parts  of  complex  numbers  z,-  of  unit 
modulus  as  above.  The  ou-diagonai  elements  of  A will  now  vanish  for  complex 
numbers  z^.  subject  to  the  same  restrictions  (with  n replaced  by  r ) as  used  for 
the  first  series.  Moreover,  since  it  has  been  shown  that  the  main  diagonal  ele- 
ments of  A must  equal  n/3  for  an  absolute  maximum,  v must  be  chosen  so  that 
ry2  + s = n/3  or  y = \f(n73^sVr.  Now  2;-i(xji)2  = r/2  = E^.jfxj;)2  so  that  the 
other  main  diagonal  elements  of  A have  the  value 


r + s n n 

~~  ~~2  6 ~ 3 

Discussion 

If  the  locations  are  to  be  usable  on  a ballistics  test  range,  the  second  rule 
for  deriving  further  optimum  locations  must  be  used  with  discretion.  Changing 
the  sense  of  a station-missile  ray  evidently  puts  the  station  initially  on  the 
upper  nappe  of  the  cone.  It  is  possible  for  several  such  changes  to  be  made  in 
such  a way  that  the  result  cannot  be  rotated  so  as  to  have  all  stations  in  the 
range  plane  and  the  missile  above  it. 

If  a configuration  has  been  modified  by  this  rule  and  successfully  rotated 
to  place  the  stations  in  the  range  plane,  the  stations  will  always  lie  on  a hyper- 
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bola  in  the  range  plane,  with  some  stations  on  each  branch.  In  a general  way, 
the  final  result  may  be  expected  to  have  some  stations  near  the  range  line  and 
some  quite  far  off-range,  so  that  locations  obtained  in  this  way  are  considered 
less  useful  for  a ballistics  test  range.  In  fact,  practical  consideration  need 
hardly  be  given  to  locations  othe"  than  those  of  the  first  series. 


THE  DOVAP  CASE 


Summary 

As  in  the  radar  case,  optimum  locations  in  the  dovap  case  may  be  specified 
in  terms  of  a circular  cone  with  the  missile  at  the  vertex  (Fig.  1).  The  distinc- 
tive transmitter-missile  ray  in  dovap  is  the  axis  of  the  cone,  for  which  the 
element-axis  angle  is  arc  cos  1/3  ~ 7CP  32'.  The  receiver-missile  rays  ore 
elements  of  the  cone  and  are  distributed  about  it  in  the  same  way  as  in  the 
radar  case:  (1)  equally  spaced  around  the  cone  or  (2)  for  n > 6,  divided  into 
groups  of  not  less  than  three  receivers,  with  the  members  of  each  group  equally 
spaced  and  the  groups  arbitrarily  oriented  with  respect  to  each  other.  Other 
spacings  indicated  in  Appendix  B apply  to  dovap  also. 

In  addition,  the  conclusion  can  be  reached  in  the  dovap  case  that  all  opti- 
mum locations  are  members  of  the  above  system.  This  is  shown  below  following 
the  nroof  that  the  locations  stated  are  optimum. 


Derivation 

For  the  dovap  case,  it  is  convenient  to  suppose  the  matrix  X to  have  n + 1 
columns.  Each  of  the  first  n columns  is  a set  of  direction  cosines  of  a receiver- 
lines  ire  ray,  cut u , .,*j  uto t Cultiutu  the  uct  of  the  direction  ccomcc  of  the 
transmitter-missile  ray.  From  the  discussion  of  the  dovap  case  in  Ref.  1,  it  is 
evident  that  the  quantities  b{/  (defined  on  page  9 of  this  reference)  are  sums  of 
direction  cosines.  Thus,  the  matrix  (6<;)  car;  be  obtained  by  multiplying  X by 
the  matrix  P on  the  right  where 


0 

0 ... 

°\ 

/ 0 

1 

0 ••• 

°\ 

: 

0 

1 ... 

0 

\° 

0 

0 • • • 

:/ 

\i 

1 

1 • • • 

1/ 

and  where  P has  n + 1 rows  and  n columns.  Then  the  determinant  A = |(6i;-)(  (>;,)'! 
of  Ref.  1 becomes  \XPP'X'\  - \XAX'\,  where 
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1 = 


Z1 

0 

0 • * • 

0 

1 

/ 0 

1 

0 . * • 

0 

1 

■ 

0 

1 ... 

0 

1 

\° 

0 

0 ... 

1 

1 

\i 

1 

1 ... 

1 

n 

\ 

/ 


has  n + 1 rows  and  columns,  and  is  symmetric  and  nonnegative.  If  *a  denotes 
the  CXth  row  of  A',  then  the  clement  in  the  ath  row  and  /3th  column  of  A = |AdA'| 
is 

t n 

xoAr  & ~ ^ xCx.ix/3i  + 2 ^xOUxfi,n  + \ + xa.,n+  1 x[id  + nx<X'Ti  + lxfl,n  + l 

i-1  «- 5 

By  Theoicm  1,  the  discussion  may  be  concerned,  without  Iosb  of  generality, 
with  the  case  in  which  the  transmitter-missile  ray  is  vertical:  (*i(„  + i,  *2,n+l> 
= (0,  1,  0).  In  this  case, 


yxA%\  = Si  = 

x2Ax2  •=  s_i  = ^ n + *2,)2 

*3/1*3  = S3  = £x& 

x\Ax2  = x2Ax\  = £ (1  + *2«)*i< 

*1-4*3  = *3  Axx  = 2*li*3« 

r-.d*3  = rs4  = V(1  + *2;)*3i 

where  the  sums,  here  and  below,  are  over  the  range  i = 1,  • • • , n. 

The  proof  that  the  locations  given  in  the  summary  of  this  section  are  opti- 
mum proceeds  somewhat  indirectly.  It  beginB  by  finding  the  maximum  of  the 
product  S1S2S3  of  the  main  diagonal  elements  of  A,  which  is  never  exceeded  by 
A itself  according  to  Theorem  2.  Since  *?,  + rf,  •*-  rf;  - 1 ( this  product  may  be 
written 

SiS2S3  = Q = (£  *u)[Z  (1  + *2«^2] [Z(1  “ *i (i  ~ *2i)] 

whence 


dO 

— = 2*hS2S3  — 2*i/;SiS2 

°x\k 

dO 

— 2(1  x2k'S\$$  — 2*2 ^'^'2 

dx2k 

Since,  obviously,  each  Sa>  0 when  (J  is  maximum,  setting  <?t//‘?*it  =0  implies 
S3  = Si,  and  setting  do/ dx^  = 0 implies  *24(^3  — S21  = — S3.  Now,  if  S3  =S2,this 
shows  S3  = 0;  hence,  for  a maximum 

(7)  *2*  = S3/(S2  “ ->3^ 

fot  each  k = 1,  2,  • • • , n.  This  means  that  all  the  *2*  are  equal,  say  to  some 
vnriable  y (-1  <y<l),  and  all  the  receiver-missile  rays  lie  in  a right  circular 
cone. 
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The  equality  of  all  the  x^kt  which  has  to  be  satisfied  in  order  for  Q to  be 
maximum,  now  simplifies  the  conditions  of  Theorem  2 which  result  in  A «»  Q. 
Thus,  the  off-diagonal  elements  are  now  zero  if 

(8)  £ *H  = Z *3i  = 0 Qnd  £ = 0 

In  addition,  it  is  shown  above  that  £,  = S3  (i.e.,  that  £ r’  = £ %3i).  This  condi- 
tion and  Eq.8  constitute  the  conditions  satisfied  by  the  optimum  locations  given 
in  the  radar  case.  Hence,  the  same  distributions  of  rays  around  the  cone  give 
optimum  locations  in  the  dovap  case.  Moreover,  since  x2’  + x2i  + *3;  = 1, 

£ *u  + £ *3,  - r.(l  - y2 1 

so  that 

Si  = S3  = (n/2)(l  - y2) 

and 


Q - A - S,S2S3  - l(»/2)(l  - y2)]Ml  + yl2  j[(n/2)(l  - y2)] 
= (n3/4)[(l  — y2)(l  + y)]2 


and  the  value  of  A is  a function  of  the  cone  angle  alone  when  the  receivers  are 
distributed  about  the  cone  as  specified. 

The  above  relations  also  determine  the  cone  angle  as  follows:  From  the 
relations  on  page  7, 

Si  + S2  4-  Sr  = n + 2 £ in;  + £ +■  x~;  + xf;)  - 2.-.Q  1 y) 

ao  that 


S-2  = 2n(l  + y)  - Sj  — S3  * nd  + y)2 

However,  by  Eq.7, 

_ _ /_  1 \ 

o2  -03^i 

2y 

since  S3  *=  (n/2)(l  — y2).  Equating  these  expressions  for  S2  yields  (1  - y)/ 2y  = 1 
or  y * 1/3.  Thus,  the  cone  angle  is  arc  cos  1/3,  and  the  maximum  value  of  A is 
256.-.V729. 


To  prove,  in  conclusion,  that  all  optimum  locations  are  included  in  those 
discussed  above,  suppose  a given  configuration  is,  in  fact,  optimum.  As  above, 
it  may  be  assumed  in  view  of  Theorem  1 that  the  transmitter-missile  ray  is 
vertical.  Now  the  argument  leading  to  Eq.7  shows  that  the  receiver-missile 
rays  must  lie  in  a circular  cone  and  that  S j = S3.  Further,  it  . l shown  above, 
for  conical  configurations,  that  the  absolute  maximum  of  Q is  256n3/729  and 
that  A = O when  Eq.8  are  satisfied;  hence,  tbq.H  are  satisfied  by  the  given 
configuration.  Thus,  if  n < 6,  the  given  configuration  is  one  of  the  type  specified 
in  the  summary  of  this  section;  if  n > 7,  other  configurations  are  possible  (shown 
in  Appendix  B). 


8 


NAVORD  REPORT  2099 


THE  HYPERDOP  CASE 


Summary 

Tfte  v.ork  on  the  hyperdop  case  has  produced  two  results  of  interest.  The 
first  is  that,  for  any  location  of  stations  for  which  all  station-missile  rays  are 
elements  of  a circular  cone,  the  methods  of  data  reduction  of  Ref.  1 fail.  The 
second  is  a set  of  sufficient  conditions  for  an  optimum  which  can  be  satisfied 
for  at  least  some  n and  for  which  the  stations  cannot  be  placed  in  the  range 
plane  with  the  missile  above  it. 

Darivatian 

The  first  result  follows,  in  fact,  directly  from  the  treatment  of  this  case 
in  Ref.  1.  To  show  this,  note  that,  by  Theorem  1,  this  cone  may  be  rotated  so 
that  its  axis  is  vertical,  resulting  in  the  equality  of  all  vertical  components, 
*2 However,  the  x 2l-  are  the  quantities  (yp  — T,-)/0j  on  page  4 of  Ref.  1 so  that 
the  differences  c^2  vanish.  Then  grs  - 0 when  either  r=  2 or  s = 2 so  that 
A «=  0,  and  the  estimate  of  elevation  of  the  missile  cannot  be  found.  Continuity 
considerations  show,  of  course,  that  the  accuracy  of  the  position  determination 
is  poor  for  any  configuration  “approximately”  conical. 

The  se  ccnd  result  is  obtained  by  using  expressions  for  the  elements  of 
A.  Evidently,  the  quantities  cty*  of  Ref.  1 are  the  differences  xk(  - xkj  so  that 
the  element  grs  in  the  nth  row  and  sth  column  of  A is 

fcrs  ^ ^ yxri  “ xrjnxsi  ~ xsj ' 

i<i 

= Z ixrixsi  + xnxsj)  - Z (xrixsj  + xnxsi' 

i< / i<j 

= (jl  — ^ xrixsi  ~ ^ xrixsj 

i 

= 'lZIr>I»i‘( Z V'i) 

with  Py  = 1 and  all  pairs  ft,  /)  with  i < / used.  It  is  easily  verified  that  the 
determinant  |grj|  is  of  the  form  \XAX'\,  where  A is  the  symmetric,  nonnegative 
matrix  whose  main  diagonal  elements  are  all  n — 1,  and  whose  off-diagonal 
elements  are  — 1.  Thus,  the  product  of  main  diagonal  elements  of  A is 

Q = [nZ  4 - (Z  *u)%  I 4 - (I  x2i)%  1 4 ~ (Z*3i)2] 

it  is  necessary  to  show  that  the  partial  derivatives  of  Q vanish  at  a maxi- 
mum. This  is  done  by  use  of  the  auxiliary  function 

= [n  Z -l2.  - (Z  *i.)2]["  Z x2i  - (Z  *;.)2J[»  Z (1  - *,2  - 4>] 

= ii  j 

in  which  no  restriction  is  placed  on  the  2 n variables  rj;,  x2i*  From  comparison 
of  the  expressions  for  Q and  O*,  it  is  evident  that,  when  4 + 4 < It  f°r 

every  choice  of  signs  for  *3,-=  ±\/\  - 4 - 4't  w'lh  O*  = Q when  X x3l-  = 0.  Now, 
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Q*  is  defined  and  has.  partial  derivatives  for  all  x2l-;  0*  < 0 for  some  x2l-; 
and  Q*  < 0 if  any  |*i,|  or  |x2t-|  *8  sufficiently  large.  Hence,  <)*  has  a maximum  at 
a finite  point,  and 

<50*/&ci*  = d<J*/dx2k  = 0 

at  this  maximum  for  all  k.  However, 

dQ*/dxlk  = 52[S3(2nxu  - 2£*lf)-  2nSi*u] 

which  vanishes  when  S2  = 0 or  when  n(S3  — Sj  )*u  = 2 xy.  Now,  if  S3  4 S],  this 
last  equality  can  only  be  satisfied  if  all  xj*  are  equal,  which  results  in  A = 0, 
as  is  shown  in  the  first  paragraph  of  the  derivation  of  this  section.  On  the  other 
hand,  if  S3  =Sj,  then  2 xj,-  = 0,  and  a similar  argument  relative  to  dO*/dx2k 
shows  that,  for  a maximum,  2 x2i-  = 0 also.  At  a maximum,  then,  Q*  reduces  to 

"3(Z  *u)(Z  X2i)(»  - 1 *?<  - I 4) 

From  the  discussion  cf  the  radar  case,  the  absolute  maximum  of  this  quantity 
is  n3(/i/3)3,  attained  whenever  2 Xj,-  = 2 x2i  = n/3. 

Thus,  if  a matrix  X exists  for  which 

2>a,  =0,  2*ai«»/ 3 (a  -1,2,  3) 

i 

and 

£ *a ix/3i  = 0 foe  ^ /3> 

i 

then,  using  Theorem  2,  t\  = Q * Q*,  and  the  absolute  maximum  of  all  of  these  is 
attained.  The  first  three  of  these  conditions  show  that  some  of  the  elements  of 
each  row  of  X are  positive  and  some  are  negative,  since  A = 0 if  all  elements 
of  any  row  are  zero;  and  it  is  easily  verified  that  the  same  statement  will  hold 
after  any  rigid  rotation.  This  means,  of  course,  that  some  stations  are  above 
the  missile  and  others  are  below  it  however  the  system  is  rotated. 

Discussion 

It  is  not  known  whether  the  sufficient  conditions  for  an  optimum  of  the 
previous  discussion  can  be  satisfied  for  all  n > 4.  It  is  easily  verified,  how- 
ever, that  configurations  having  the  missile  at  the  center  and  the  stations  along 
the  rays  to  the  vertices  of  a regular  tetrahedron  and  a cube,  respectively,  satisf) 
these  conditions.  Obviously,  configurations  built  up  from  tetrahedra  and  cubes 
in  arbitrary  relative  orientation  are  optimum  also. 


NUMERICAL  ILLUSTRATION 

A number  of  numerical  evaluations  of  the  determinant  were  made  for  the 
radar  and  dovap  cases  in  order  to  study  their  variations  as  the  missile  approached 
the  optimum  and  as  the  number  and  arrangement  of  stations  were  changed.  Tht 
same  missile  positions  were  used  for  both  radar  and  dovap,  including  four  points 
on  a parabola  passing  through  the  launcher  with  vertex  50,000  feet  down-ra  ige 
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and  20,000  feet  in  elevation.  Except  for  one  case,  the  configurations  are  opti- 
mum for  the  missile  at  the  vertex,  and  they  differ  between  the  radar  and  dovap 
cases  only  in  the  necessarily  larger  circle  in  the  range  plane  for  the  dovap 
receivers. 

The  coordinates  of  the  missile  position  are  given  below,  using  x for  down- 
range,  y for  elevation,  and  z for  off-range. 

Missile  position,  ft 

2 3 4 

45,000  30,000  5,000 
19,800  16,800  3,800 
0 0 0 

Five  arrangements  of  stations  were  used.  The  first  had  the  minimum  number 
(3)  of  stations;  the  second  and  third  had  4 stations  in  different  orientations  with 
respect  to  the  down-range  axis;  the  fourth  had  6 stations  at  the  vertices  of  a 
regular  hexagon;  and  the  fifth  had  3 stations  in  a nonoptimum  configuration 
simulating  a case  of  4 stations  with  no  record  from  one  of  the- in.  The  coordinates 
of  these  configurations  are  given  in  Table  1. 

TABLE  1,  Configurations  of  Stations  Used  in  Numerical  Calculations 


The  dovap  transmitter  always  has  coordinates  (50,000,  0,  0),  and  y is  always  zero. 


Configuration 

Coordinate,  ft 

Radai 

Dovap 

X 

2 

X 

2 

1.  Equilateral  triangle 

78,284 

0 

106,569 

0 

35,858 

24,495 

21,716 

48,990 

35,858 

-24,495 

21,716 

-48,990 

2.  Square  with  aiag.  on  range  line.... 

•*«  nn  4 

« i 

0 

106.569 

0 

50,000 

28,284 

50,000 

56,569 

21,716 

0 

-6,569 

0 

50,000 

-28,284 

50,000 

-56,569 

3.  Square  with  2 sides  parallel 

70,000 

20,000 

90,000 

40,000 

to  range  line 

30,000 

20,000 

10,000 

40,000 

30,000 

-20,000 

10,000 

-40,000 

70,000 

-20,000 

90,000 

-40,000 

4.  Hexagon  

73,284 

0 

106,569 

0 

64,142 

24,495 

78,284 

48,990 

35,858 

24,495 

21,716 

48,990 

21,716 

0 

-6,568 

0 

35,858 

-24,495 

21,716 

-48,990 

64,142 

-24,495 

78,284 

-48,990 

5.  Isosceles  right  triangle0. 

78,284 

0 

106,569 

0 

50,000 

28,284 

50,000 

56,569 

21,716 

0 

-6,569 

0 

° Simulating  configuration  2 with  one  station  missing. 


Coordinate  l 

x 50,000 

y 20,000 

z 0 
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TABLE  2.  Values  of  Determinants  Obtained  in  Calculations 


Station 

configu- 

ration 

Radar  missile  position 

Dovap  missile  position 

1 

2 

3 

4 

1 

1 

2 

3 

4 

1 

1.000 

0.975 

0.555 

0.005 

9.481 

9.050 

4.677 

0.084 

2 

2.370 

2.353 

2.126 

0.033 

22.475 

21.6C4 

13.262 

1.752 

3 

2.370 

2.352 

1.767 

0.018 

22.475 

21.601 

12.731 

0.366 

4 

8.000 

7.881 

6.968 

0.104 

75.852 

72.913 

44.417 

5.480 

5 

0.592 

0.588 

O.I.ji 

0.008 

5.619 

5.389 

3.201 

0.376 

The  determinants  were  evaluated  for  every  combination  of  missile  position 
and  station  configuration.  T*>>*  values  found  are  ^iven  in  Table  2.  The  following 
remarks  may  be  made  relative  to  the  values  in  this  table: 

1.  The  values  do  not  permit  a general  comparison  of  accuracy  between 
radar  and  dovap,  since  account  is  not  taken  of  other  sources  of  error. 

2.  Dovap  consistently  loses  more  accuracy  than  radar  in  missile  positions  2 
and  3 but  less  in  position  4. 

3.  Of  the  two  4-station  configurations  (2  and  3).  the  one  with  2 stations  on 
the  range  line  is  better  for  both  radar  and  dovap. 


FURTHER  OBSERVATIONS 

The  determinants  A for  optimum  locations  on  cones  with  vertical  axes  have 
properties  worthy  of  further  discussion.  The  fact  that  the  off-diagonal  elements 
are  always  zero  mc.ins  that  the  three  components  of  the  fix  are  uncorrelated. 
This  property  simplifies  further  statistical  inferences  based  on  the  data,  such 
as  estimates  of  velocity  and  acceleration.  Again,  the  main  diagonal  elements  of 
A remain  in  fixed  ratio  for  all  n in  each  type  of  instrumentation.  These  ratios 
show  that  (1)  for  radar,  the  optimum  ellipsoid  of  concentration  is  always  a sphere 
and  (2)  for  dovap,  the  vertical  axis  of  the  ellipsoid  is  always  one-half  the  other 
axes,  which  are  always  equal. 

Another  criterion  of  a :curacy  having  strong  intuitive  appeal  is  the  expected 
square  of  the  radial  error  of  the  fix.  The  results  of  this  report  suggest  that,  for 
radar,  the  same  locations  are  optimum  hv  this  criterion.  For  dovap,  however, 
this  criterion  apparently  leads  to  different  locations  for  which  the  components 
of  the  fix  may  not  be  ur.correlated;  nor  is  there  good  reason  to  suppose  that  the 
ellipsoid  of  concentration  is  spherical.  Locations  optimum  by  the  alternative 
criterion  should  not  be  markedly  different  from  those  of  this  report,  however, 
in  the  sense  that  regions  of  “good”  fixes  by  the  two  criteria  should  overlap 
extensively. 
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Appendix  A 

PROOFS  OF  THEOREMS 

The  twr  theorems  given  in  the  Introduction  and  a third  theorem  of  general 
character  are  proved  in  this  appendix.  Theorem  3 actually  led  to  the  discovery 
of  optimum  locations  in  the  dovap  case,  although  the  results  are  more  easily 
given  otherwise. 

Theorem  1.  Rigid  rotation i in  space  of  the  system  of  rays  whose  direction 
cosines  constitute  X leave  A = unchanged. 

PROOF.  Let  R be  the  matrix  oi  the  rotation  so  that  RX  is  the  matrix 
of  direction  cosines  after  the  rotation.  Then  R is  orthogonal,  \R\  **  1,  and 

\RXAX'R'\  - \R\  \XAX'\  |/T|  = \XAX'\. 

The  proof  of  Theorem  2 requires  a theorem  of  Hadamard.  This  is  given  on 
pages  219  and  220  of  Kef.  4 as  follows:  Let  D = | | be  an  n-row  determinant 
of  real  elements,  and  let  X".j  .7*  = p,  > 0.  7 her.  the  maximum  of  D2  is  pj  • • • p„, 
and  this  maximum  is  attained  if  and  only  if 

Z = 0 

k-i 

UtEOREM  2.  Let  X = (x:j)  he  a matrix  of  variables  with  m rows  and  n columns 
of  rank  m,  and  let  *i  *=  t*n,  * * * , xin ) be  the  ith-row  vector  of  X.  If  the  matrix  A 
is  symmetric  and  positive  semidefinite  (i.e.,  all  its  orincipal  minors  are  non * 
negative),  aid  if  x,/4xj  > 0 for  i = 1,  • • • , m,  then 

\XAX'\  <{XlAX\)  • • • (xmAx'J 

where  the  equality  sign  holds  if  and  only  if  AxJ  = 0 when  i ^ /. 

PROOF.  Let  B = ( 6,-j- ) be  a square  matrix  such  that  /?/?'  = A,  and  set 
X3  = i ’ - (yij).  If  m < n,  let 


where  the  Zjh  are  chosen  as  solutions  of  the  homogeneous  equations 
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Z y i'-Sjk  = 0 

*-i 


(i  = 1,  • • • , m;  /'  = m + 1,  • • • , n) 


which  are  orthogonal  and  normalized — that  these  exist  may  be  seen  from  the 
theory  of  linear  spaces — so  that 


Z */*-!.  Z *»**/»  "° 

*-i  k- 1 


(i  4 j;  i,  j > m + 1) 


and  let  Y = ) il  m = n.  Now 
On  the  other  hand. 


YY ' 0 

0 


= | YY' 


|T|2  <(Z  yu)  • • • (Z  ym*) 

by  the  Iladamard  theorem,  and  the  equality  holds  if  and  only  if  X y,*y;*  =0  for 
i 4 /.  Generally, 

Z ‘jikYjk  = Z (Z  xisbskj( Z 
“ Z (Z  bskbtk)*uxit 

= Z astxisxjt  = 

Thus,  for  f «*  /, 

|yr|  = |^r|  < . . . (*„,/!*;) 

while  the  equality  is  seen  to  hold,  by  the  Hadamard  theorem,  if  and  only  if 
XjAx'j  = 0 for  t ^ /.  This  completes  the  proof. 

Note.  \XAX'\  = i > 0. 

Theorem  2,  when  applied  to  A = shows  that 

A ^ = (Z  ad*l.*l/)(Z  aHx2ix2, “.>*3i*3/) 

while  A = Q (4  0)  if  and  only  if  = 0 for  (X  4 fi~ 

Theorem  3,  which  is  given  below,  simplifies  the  results  obtained  by  applying 
the  method  of  Lagrange  to  the  problem  of  maximizing  Q under  the  restrictions 

<91  *1.  + x2i  + x3i  - 1 

lue  method  yields,  in  addition  to  Lq.9,  the  necessary  conditions 


00) 


Z + ^k*lk  - 0 

; 

Z a*l*2#  + ^2*  = 0 

j 

Zav;*3,  + = 0 
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for  k = 1,  • • • , n,  where 

f = W - VS1S3.  a<a3>  = a*/s,s2 


(11) 


ai)x0Lix 


ix0Lj 


(a  - 1,  2,  3) 


THLORfcM  3.  Necessary  and  sufficient  conditions  that  a matrix-  A constitute 
a solution  of  the  systemlO  subject  to  Eq.  11  are  that  there  exist  a matrix  A = diag 
(Aj,  • • • , A„)  and  constants  p,-  > 0,  i = 1,  2,  3,  such  that 

(12)  (A  + PiMx(  = 0 

and 


(13)  XiNx'i  = -(1/VPiP2P3^ 

for  i = 1,  2,  3,  where  Xj  is  tAc  row  vector  (*n*  * * * * *<«)• 

PROOF.  For  the  proof  of  the  necessity  of  the  conditions,  if  A is  a solution 
of  system  10,  the  matrix  A and  the  constants  p,  exist  such  that  Eq.  12  is  satis- 
fied. Moreover, 


Pi  = I/S2S3*  P‘i  ~ 1/S,S3,  p 3 = 1/SjS2 
From  these  relations,  it  follows  easily  that 

(14)  i'l/pi  ••  S2/p2  “ S3/p3 

It  also  follows  from  Eq.  12  that  x,(A  + PjA)*j  *»  0;  that  is, 

(15)  S,  - = -piXiNx'i 

and  substitution  in  Eq.  14  resuPs  in 

(16)  xjAx'i  «=  x2Ax2  =••  x3Ax3 
Now  the  corresponding  value  of  Q,  say  Qq,  is  given  by 

On  = (— n, r , A vIVS-.S- 

whence,  since  PiA'2.S3  = 1, 

(17)  Qn  =—  x\hx\ 

Out 


$1  /^2\^3  ^0 

Pi  \P2/  P3  P1P2P3 

= (— X|Axi)(— x2Ax2)(--x3Ax3) 

= ^03 

that  is, 

(18)  (J0  = I/VP1P2P3 

and,  using  Eq.  16  and  17,  it  is  seen  that  Eq  . 13  is  established. 

For  the  proof  of  the  sufficiency  of  the  conditions,  it  follows  from  Eq.  13  that 
(x2Ax2)(x3Ax3)  = l/p,p2p3 

whence 


Pi  A ^ A/(-p2x2Ax2)(-p3x3Ax3) 
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0y  Eq.15,  this  relation  is  equivalent  to 

Pi  A = A/S2S3 

and  similar  expressions  hold  for  pjA  an(i  P3A.  Thus,  with  = PoXki  Eq.l2and 
13  imply  that  X satisfies  a system  like  the  system  10  subject  to  Eq.ll,  and  the 
proof  is  complete. 

Note.  It  follows  from  Eq.  9,  16,  17,  and  18  that 

3<?0  “ ■M*?l  + *21  + *3 1‘) 

“ “Z  *1  ” 3/v/ri?2P3 
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Appendix  B 

METHODS  OF  CONSTRUCTING  FURTHER 
CONFIGURATIONS  FOR  LARGE  n 

It  is  shown  in  the  body  of  the  report  that  necessary  and  sufficient  conditions 
for  a set  of  plane  direction  cosines  *3*)  to  determine  an  optimum  location 

o it koe  rsdur  or  dovor'  ors 

Y *ik  = n/2  Z*i**3*=°  (4=1,  •••,») 

It  is  also  shown  that  the  monic  complex  polynomial  whose  zeros  are 

z*  - -t-  *3*^-!.  4 = 1,  • • • , n, 

is  of  the  form 

(20)  zn  + a3Zn~3  + • • • + an_jz  + on 

Conversely,  it  is  clear  that,  for  any  complex  polynomial  of  the  form  20  with 
zeros  of  modulus  unity,  an  optimum  location  can  be  constructed.  As  a con- 
sequence of  the  following  theorem,  the  polynomial  is  seen  to  take  a still  more 
special  form. 

TltUJREM  4.  if  the  complex  polynomial 

zn  r axzn~l  + a1zn~'i  + ••»  + «„ 

with  zeros  zj,  • • • , zn  is  such  that  |z*  | = ).  for  all  k = 1,  • • • , n,  then  aj  = ano„_,-. 
PROOF.  As  is  well  known, 

ai  = (-l)'  Y zk\  * * * zkj 

fc, < 

~ (-1V  2 • • • zkj) 

*!<•••<*/ 

since  |z*|  = 1.  Hence, 

afan  ■=  (-l)/+n  2 (^1*2  • ’ ’ zr,)/(*k  1 * * • zkj) 

kx<‘--<kj 

The  right  member  is  now  just  an-j,  and,  since  |a„|  = 1,  the  theorem  follows. 

From  this  theorem,  it  follows  that  optimum  locations  are  determined  by  com- 
plex polynomials  of  the  form 

(21)  zr  + u3Zn-3  + • • • + a$z3  + 1 

witii  zeros  of  unit  modulus,  where  the  specialization  a„  = 1 merely  rotates  the 
zeros  equally  around  the  unit  ircle.  The  only  polynomials  of  .his  form  for 
n = 3,  4,  5 are  z"  + 1,  and  the  zeros  of  these  are  equally  spaced  on  the  unit 
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circle;  hence,  the  only  optimum  locations  of  3,  4,  or  5 receivers  have  them 
equally  spaced  around  the  cone.  For  n - 6,  the  polynomial  is  quadratic  in  z3; 
hence,  each  of  its  factors  ieads  to  3 stations  in  an  equilateral  triangle,  while 
the  triangles  may  be  arbitrarily  orie:  ed  with  respect  to  each  other. 

For  larger  values  of  n,  various  special  methods  have  been  used  to  find 
further  optimum  locations.  For  example,  using  the  above  theorem,  if  a poly- 
nomial of  the  form 

zr  + b^z'-2  + • • • + 6r-2*2  + 1 

has  zeros  of  unit  modulus,  then  the  polynomial  obtained  from  it  by  replacing 
z by  z1  for  s > 2 is  a polynomial  of  the  form  21  and  ieads  to  an  optimum  loca- 
tion for  n = rs  stations.  This  method  does  not  lead  to  further  optimum  locations 
for  r < 5 or  s > 2. 

Another  special  method  uses  Eq.  19  in  their  equivalent  trigonometric  form 
2_  cos  0*  = 0 sin  W*  = 0 

£ cos  2 0*  = n/2  2^  sin  0*  cos  0*  = 0 

where  0*  is  the  argument  of  the  complex  number  z*.  This  method  has  led  to 
additional  configurations  for  n = 7 in  the  special  case  04  = -0i,0j  = -0 2> 
06  = -0- , and  07  « 0.  This  assumption  reduces  the  above  equations  to 

f 1 4-  f2  + *3  — —■  ^ 

,2  ^ ,2  . ,2  5 

<1  + t2  + *3  = 4 

where  **  = cos  0A,  4 = 1,2,  3.  These  equations  have  the  solutions 
*2  = - \ (2*i  + 1 - V9-4*,  -12*2) 
t3  =-1(2*,  + l+V9-4*,-12*2) 

ft Oiii  which  aii  locations  may  be  obtained  by  choosing  *i  from  the  interval 
[i,  (2x/7  - l)/6],  L ocations  obtained  in  this  way  include  the  following: 

(*1  «•  y).  0j  = 60°,  02  = 9CP,  03  m 180°;  an  equilateral  triangle  with  a vertex 
at  —1  and  a square  with  diagonal  on  the  real  axis 

*1  = cos  (360°/7).  A regular  heptagon 

*1  = ^\/5.  0i  = 45°,  02  = 120'’,  03  = 135°;  another  case  of  equilateral  triangle 
and  square 

*1  = (2\/T  — l)/6.  0i  = 44°  20'  02  - 03  = 127°  25';  an  irregular  pentagon  with 
2 stations  at  each  of  2 vertices 
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NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  a DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS  , 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OP-  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTUR1  . 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANT  WAY  BE  RELATED  THERETO. 


WAY  BE  RELATED  THERETO. 
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